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*~^' The article by Hao and Scheraga [1] states that it is the first application of the so-called 

"entropic sampling Monte Carlo (ESMC)" method [2] to the protein folding problem: "It 
r-j \ should be noted that the earlier procedures are all aimed at simple systems such as spin 

^ ' lattices and fluids" . However, there already exists an earlier work that introduced the same 

method to the protein folding problem [3]. In the latter article the method was referred 
rS\ to as "multicanonical Monte Carlo (MMC)" [4-7]. Although MMC was quoted in Ref. [1], 

Hao and Scheraga failed in recognizing the identity of ESMC to the earlier MMC method. 
The reason is possibly that the authors of Refs. [1,2] seem not to have been aware of 
Refs. [6,7], where the idea and numerical methodology had been more explicitly stated 
than in Refs. [4,5]. These latter works (Refs. [6,7]) had already appeared in print before 
Ref. [2] was submitted for publication. The erratum of Ref. [2] corrects for this neglect, but 
Hao and Scheraga [1] do not refer to this erratum. It is the purpose of this comment to 
give a simple proof that the two superflcially differently looking approaches are identical. 
Essentially, we follow some earlier private correspondence of one of the authors with Lee 
[8]. 

In ESMC configurations with energy E are updated with a weight 

w{E) oc e-^(^) , (1) 
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(see, for example, Eq. (3) of Ref. [1]) where S{E) is the microcanonical entropy. We have 
set the Boltzmann constant k equal to 1 for simplicity. On the other hand, in MMC a 
simulation is performed in a multicanonical ensemble [5,7], which in turn is defined hj the 
condition that the probability distribution of the energy shall be constant: 

P{E) = n{E)w{E) = const , (2) 

(see, for example, Eq. (4) of Ref. [6]) where n{E) is the spectral density (or density of 
states) . Hence, this condition implies a weight factor 

w{E)<xn-\E) . (3) 

It is standard textbook knowledge of statistical mechanics that the density of states n{E) 
and the entropy S{E) are related by 

n{E) = e^(-^) . (4) 

Hence, ESMC and MMC have the identical weight factor w{E) (see Eqs. (1), (3), and 
(4)). Both methods are based on the same idea that a uniform (flat) distribution in energy 
is obtained (rather than the usual bell-shaped canonical distribution) so that even regions 
with small n{E) or S{E) may be explored in detail and thus a good estimate of the density 
of states or entropy may be obtained [2,4-7]. Therefore, the two methods are conceptually 
identical. 

The identity of the two methods was obscured possibly because most works based on 
MMC use the following parametrization of the weight factor of Eq. (3) : 

«;(£;) oce-^(^)^-«(^), (5) 

which is different from that of ESMC (see Eq. (1)). Here, P{E) can be interpreted as an 
"effective" inverse temperature. While not necessary, this choice is particularly appealing 
when one wants to combine canonical and multicanonical simulations [3-7]. At flrst sight, 
one may be puzzled by the fact that ESMC has only one parameter S{E) , while MMC seems 
to have two parameters a{E) and (3{E). However, there is also only one free parameter 
for MMC, since a{E) is obtained from (3{E) by the relation [4,6] 

a{E) = a{E') + {f3{E') - (3{E))E , (6) 

where E' and E are adjacent bins in the histogram. Both ESMC and MMC discretize the 
energy space with a certain bin size and consider the (unnormalized) histogram H{E) of 
the energy distribution. The parameters S{E) for ESMC and (3{E) for MMC are both 



determined by an iteration of simulations with small Monte Carlo statistics. In ESMC 
S{E) is updated by [2] 

S'+\E) = S'{E) + \nH\E) for H\E) > 1, (7a) 

and 

S'+\E) = S'{E) for H'{E) = 0. {7b) 

Here, i is the iteration number and H^{E) is the histogram of the energy distribution after 
i-th iteration. 

We now show that the recursion relation, Eq. (7a), for ESMC is identical with that 
for MMC derived earlier by Berg and Celik [6] . The identity of the weight factors of both 
methods implies the following relation among the parameters at the i-th iteration (see Eqs. 
(1) and (5)): 

S\E)=f3'{E)E + a\E). (8) 

By Eq. (6) we can eliminate a^{E) and obtain 

S\E') - S\E) = f3\E'){E' - E) . (9) 

This means that once we have ESMC parameter S{E), we can obtain MMC parameters 
l3{E) and a{E), and vice versa. Indeed, Eq. (9) with the recursion relation Eq. (7a) gives 
the following relation 

/J«(E')^/3'(E0 + ^..(^), (10) 

which is exactly the recursion relation given in Eq. (5) of Ref. [6] {E' — E = —4 there). 
Therefore, the two methods are algebraically and numerically identical. Step by step the 
computer will generate the same configurations, provided that trivial details like random 
numbers are also chosen to be identical. This is why we continue to refer to the method 
as MMC in our subsequent papers [9-11]. 
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